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A  WEIGHTED  LINE-FINDING  ALGORITHM 


I.  INTRODUCTION 


This  paper  presents  two  new  algorithms  to  detect  a  line  in  a  digitized  picture.  We 
restrict  ourselves  here  to  the  simplest  case,  a  digitized  picture  with  binary  black  and  white 
pixels,  although  detecting  colinear  points  in  a  general  picture  is  possible  with  some 
modifications  in  the  algorithms. 

The  ultimate  use  of  such  algorithms  could  be: 

1.  To  take  digitized  imagery  into  a  linearized  reduced-data  format  that 
could  be  rapidly  compared  to  an  existing  digital  map 

2.  To  speed  up  automatic  line-following  procedures  that  aid  photo  inter¬ 
preters,  and  in  general 

3.  To  facilitate  the  process  of  using  image  data  in  mapping. 


In  [  1  ]  Duda  and  Hart  described  a  method  due  to  Hough  (  “Hough  transform”  ) 
that  replaces  finding  colinear  points  by  an  equivalent  problem  of  finding  intersections 
in  the  sinusoidal  curves  they  generate  by  that  transform.  That  is,  each  pixel  (  Xj,  Ys  ) 
in  the  image  space  is  transformed  into  an  r -9  parameter  space  curve: 

(1)  r  =  Xj  cos(0)  +  Y;  sin(0)  -R  =  <r<=R  0  =  <6<r 

The  parameter  space  curves  that  correspond  to  three  or  more  colinear  points,  pass 
through  a  single  point  (r^  0j),  which  corresponds  to  the  line  that  connects  the  points. 


Section  2  presents  the  new  line-finding  algorithm  (LF).  In  LF  taking  each  pair 
of  image  points  defines  a  line  and  this  line  is  transformed  into  a  point  (rj,  0j)  in  the 
parameter  space.  In  this  the  algorithm  differs  from  the  Hough  transform  algorithm  which 
operates  on  individual  image  points  taking  them  into  parameter  points.  The  new  al- 
goritnm  has  the  major  advantage  of  detecting  perturbed  lines  (lines  which  are  broken 
into  several  parts  or  have  portions  shifted  slightly  from  the  original  position)  as  entities. 


Section  3  presents  an  extension  of  the  line-finding  algorithm,  called  the  weighted 
line-finding  algorithm  (WLF).  By  using  WLF  it  is  possible  to  define  a  specific  criterion 
for  a  searched  line,  for  example  distinguishing  connected  lines  from  disconnected 
(dashed)  ones,  or  locating  a  feature  involving  a  jagged  line. 


II.  LINE-FINDING  ALGORITHM 


In  this  algorithm  each  pair  of  image  points  defines  a  line  and  this  line  is  trans¬ 
formed  into  a  point  (rj ,  0j)  in  the  parameter  space.  Here  r;  is  the  distance  to  the  origin 
from  the  line,  and  0t  is  the  angle  between  the  normal  from  the  origin  to  the  image  line, 
and  the  x-axis. 


1.  Algorithm 


Consider  a  line  passing  through  a  pair  of  points  (Xj ,  Yj )  and  (X2,  Y2).  It  is  shown 
in  Appendix  A  that  the  distance  of  the  origin  from  the  line  will  be 


(2)  I  r  1 


|X1-Y2-X2-Y,1 
V(X,  -  X2)2  +  (Yj  -  Y2)2 


The  main  steps  in  the  analysis  based  on  fig.  2,  are  that  the  distance  between  the  origin 
and  the  intersection  of  the  line  and  the  x  -axis  is 


(3)  OA  =  (X2Y,  -  XjY2)  /  (Y,  -  Y2) 
and  the  resultant  transformation  parameters  are: 

(4)  6  =  f  +  arctan  [  (  Y2  -  Y2  )  /  (  X,  -  X2  )  ) 

+  |  r  |  ifOA  <  0  or  6  <  f 

(5)  r  = 

-  |  r  |  otherwise 


To  implement  this  algorithm,  the  parameter  space  is  quantized  into  cells  and  an 
accumulator  is  assigned  to  each  cell  in  the  two-dimensional  array.  This  quantization 
can  be  confined  to  the  region  0  <  d  <  rr,  -R  <  r  <  R.  Points  outside  this 
parameter-space  rectangle  correspond  to  lines  in  the  image  that  do  not  cross  the  retina. 
R  is  the  radius  of  the  retina,  the  smallest  circle  centered  at  the  origin  surrounding  the 
image  (see  fig.  1 ). 
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Each  pair  of  points  in  image  [  (X,Y)  J  space  creates  a  corresponding  point  in  para¬ 
meter  space  counted  by  a  corresponding  accumulator  increment.  After  processing  all 
image  pixels,  the  program  inspects  the  array  to  locate  cells  with  high  counts.  For  picture 
region  edges  that  digitized  into  straight  lines  on  the  scanning  retina,  the  significant  of 
a  higher  cell  count  is  a  larger  number  of  colinear  image  points.  The  relation  between  a 
cell’s  count  and  the  number  of  colinear  points  in  the  image  is  not  a  linear  function.  In 
fact,  it  acts  so  that  those  edges  which  are  perturbed  (broken  into  several  parts  by  dig¬ 
itization)  are  smoothed  by  transformation  and  subsequent  inverse  transformation.  This 
is  the  main  advantage  of  LF  transformation  over  the  Hough  transformation. 


2.  Example  1 


To  compare  the  Hough  algorithm  with  LF  algorithm,  they  were  both  implemented 
by  the  first  author  on  the  PDP/1 1  computer  at  the  UCLA  Computer  Science  Department, 
and  applied  to  the  data  shown  in  fig.  3,  a  20  x  20  picture  with  many  colinear  points. 


Quantizing  6  at  k  =  20  nine-degree  increments  in  6,  and  quantizing  r  into 
1  =  28  one-unit  increments  in  r,  yields  the  two-dimensional  accumulator  array  of  Table 
1  for  the  Hough  algorithm,  and  Table  2  for  LF  algorithm. 

Thresholds  selected  for  the  two  techniques  permitted  detection  of  line  #4  in 
fig.  3,  an  unperturbed  edge.  With  the  Hough  technique,  the  count  in  the  parameter 
space  corresponding  to  line  #5  was  lower  than  that  of  line  #4,  and  line  #5  was  not 
detected,  although  it  is  relatively  a  long  line.  However,  LF  transformation  had  a 
corresponding  count  sufficiently  large  for  this  line  to  be  detected.  In  other  words,  LF 
transformation  has  the  significant  advantage  of  enhancing  relatively  long  perturbed  lines 
in  digitized  pictures. 


Let  us  now  consider  an  image  of  n  by  m  pixels,  and  define  k  different  quanti¬ 
zation  values  for  6.  In  the  Hough  algorithm,  the  accumulators  would  be  incremented 
kmn  times,  while  LF  algorithm  would  increment  the  accumulators  mn(mn-l)/2 
times.  The  time  required  to  process  the  sample  picture  in  fig.  3  by  LF  algorithm  was 
twice  that  of  the  Hough  algorithm,  for  k  =  20  and  m  =  n  =  20.  But  notice  that  in 
LF  transformation,  the  processing  time  necessary  is  independent  of  k.  This  is  a  second 
major  advantage  of  LF  algorithm,  since  it  allows  the  user  to  have  higher  resolution  for 
r  or  9  without  substantial  increase  in  processing  time. 
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III.  WEIGHTED  LINE-FINDING  ALGORITHM 


The  algorithm  introduced  in  section  2  needs  a  large  amount  of  processing  time  for 
large  m  or  n,  since  the  computation  required  is  of  order  (m  •  n)2 .  To  compensate 
this  weakness,  we  present  its  extension,  the  weighted  line-finding  (WLF)  algorithm.  Take 
a  point  (Xj,  Y.)  in  the  image  and  put  a  (2p  +  1)  by  (2p  +  1)  window  centered  there. 
By  giving  different  weights  to  the  lines  connecting  the  point  (Xj,  Yp  with  other  points 
inside  the  window  and  moving  the  window  through  the  image,  we  can  find  a  new  ver¬ 
sion  of  the  corresponding  parameter  space.  Through  this  algorithm,  we  can  define  a 
special  criterion  for  lines  under  search  by  defining  different  size  windows  and  different 
weighting  factors  for  neighboring  points.  The  dimensions  of  the  window  are  usually 
small,  this  limits  the  processing  time  of  the  algorithm. 


In  this  algorithm,  the  accumulators  are  incremented  an  order  of  mn  times,  which 
is  independent  of  k,  the  number  of  quantization  levels  for  6.  To  implement  this  al¬ 
gorithm,  for  each  point  (Xj,  Yj)  in  the  image  space,  we  will  consider  all  of  the  neighbor¬ 
ing  points  within  a  2p  +  1  by  2p  +  1  window  centered  at  (Xi5  Yj).  By  defining 
weights  Wkl  for  each  point  (Xt  +  k,  Yj  +  «,)  inside  the  window,  we  can  find  the 
parameters  of  the  line  passing  through  these.  Then  by  incrementing  the  corresponding 
accumulator  in  the  parameter  space  by  Wkl,  we  can  process  all  image  pixels.  The 
program  inspects  the  array  to  locate  cells  with  high  counts. 

To  compare  this  algorithm  with  the  Hough  algorithm,  WLF  algorithm  was  im¬ 
plemented  by  the  first  author  on  the  PDP/10  and  applied  to  data  shown  in  fig.  3.  By 
defining  a  7  x  7  window  and  assigning  a  weighting  factor  of  4  for  the  first  nearest 
neighbors  of  its  center,  a  factor  of  2  for  its  second  nearest  neighbors,  and  a  factor  of  1 
for  its  third  nearest  neighbors,  the  same  results  have  been  obtained  as  the  Hough  al¬ 
gorithm.  The  time  required  to  process  the  sample  picture  of  fig.  3,  by  the  two  algorithms, 
was  the  same,  and  remains  the  same  when  the  dimensions  of  the  image  are  increased. 
Table  3  shows  the  contents  of  accumulator  array  for  WLF  algorithm. 

Since  each  accumulator  of  parameter  space  will  be  incremented  if  and  only  if 
there  is  a  corresponding  line  in  the  image,  the  number  of  zero-value  accumulators  in 
WLF  algorithm  is  higher  than  that  of  the  Hough  algorithm.  This  means  that  we  need 
less  storage  to  store  the  output  data  if  we  use  a  data  structure  that  permits  us  to  store 
only  non-zero  output  data.  The  second  advantage  of  this  algorithm  is  the  capability 
of  the  algorithm  to  define  and  use  special  criteria  for  search  for  lines  through  intro¬ 
duction  of  the  weighting  factors.  Moreover,  the  processing  time  is  independent  of  k, 
and  is  an  order  of  mn.  It  thus  allows  the  user  to  have  higher  resolution  for  r  or  6, 
without  substantial  increase  in  incrementing  time.  For  an  m  by  n  image  and  2p  +  1 
by  2p  +  1  window,  the  number  of  accumulator  increments  is  almost  equal  to 
4p  (p  +  1)  mn. 


1 .  Example  2 

The  following  example  illustrates  some  advantages  of  WLF  over  the  Hough  al¬ 
gorithm  using  fig.  4  as  the  input  data.  In  this  image  line  #5  is  a  disturbed  line  longer  than 
line  #4,  and  line  #1  is  a  disconnected  line  longer  than  line  #2.  It  we  apply  the  three 
algorithms,  i.  e.  LF,  WLF  and  the  Hough  algorithm,  and  choose  different  threshold 
values  sufficient  to  detect  line  #2  in  each  case,  we  see  that  the  Hough  algorithm  will 
detect  line  #1,  #2,  #3,  #4,  and  line  #5  (line  #3  and  line  #4  will  have  the  same  parameters 
due  to  quantization  effect  in  parameter  space).  Using  LF,  all  the  lines  from  #1  to  #5 
were  detected.  The  disconnected  lines  are  distinguished  from  the  connected  ones  by 
WLF  where  only  lines  #1,  #3,  and  #4  were  detected.  Thg  output  results  of  the  three 
algorithms  are  shown  in  Tables  4,  5,  and  6.  The  count  for  line  #1  there  is  less  than  that 
of  line  #2,  because  of  the  discontinuity  present  in  line  #1,  this  occurs  although  line 
#1  is  longer  than  line  #2. 


With  m  =  n  =  30  for  fig.  4,  and  using  k  =  20,  the  processing  time  for  WLF  was 
almost  the  same  as  that  of  the  Hough  algorithm,  but  the  time  required  for  LF  was 
five  times  higher.  As  shown  in  Table  4,  5,  and  6,  the  number  of  non-zero  entries  in 
Table  4  is  higher  than  that  of  Table  5  and  6,  so  to  store  them,  the  former  needs  more 
memory.  Table  6  has  fewer  non-zero  entries  than  the  others.  This  shows  an  important 
advantage  of  WLF  algorithm. 


2.  Trade-off  among  LF,  WLF  and  the  Hough  algorithms 

The  three  algorithms,  LF,  WLF  and  the  Hough  transform,  were  compared  and 
the  results  have  been  tabulated  in  Table  7.  For  an  image  of  m  by  n  pixels,  the  com¬ 
putation  time  is  proportional  to  mn  and  kmn  for  WLF  and  the  Hough  algorithm 
respectively,  while  for  LF,  the  time  is  proportional  to  mri(mn-l)/2.  This  is  the  main 
weakness  of  LF  transforms. 


In  the  Hough  algorithm,  processing  time  is  proportional  to  the  number  of  quati- 
zation  levels  for  0,  while  this  is  not  the  case  for  WLF  and  LF.  This  makes  it  possible 
to  have  a  higher  resolution  in  parameter  space  without  higher  processing  time. 


The  storage  requirement  for  WLF  is  less  than  that  of  LF  and  the  Hough  trans¬ 
form  because  there  are  more  zero-value  accumulators  for  WLF  results. 


With  the  WLF  algorithm,  it  is  possible  to  define  an  arbitrary  criteria  for  search 
for  a  line.  By  this  capability,  it  is  possible  to  distinguish  dashed  lines  and  perturbed  lines 
from  connected  ones.  Although  with  the  LF  transform  it  is  not  possible  to  define  a 
special  criteria  for  search  for  lines,  perturbed  lines  are  enhanced  by  it.  With  respect  to 
connected  ones,  this  is  an  advantage  of  LF  with  respect  to  the  Hough  transform. 
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In  the  Hough  transform  and  LF,  the  accumulator  counts  will  show  the  number 
of  pixels  in  each  line,  while  in  WLF,  it  is  possible  to  associate  the  length  of  the  line  to 
the  count  of  accumulators.  This  makes  it  possible  to  eliminate  the  sensitivity  of  the 
transform  to  the  directional  orientation  of  the  lines. 


IV.  CONCLUSION 


In  section  2,  we  discussed  a  new  technique  to  detect  edges  in  an  image,  called  LF, 
that  is  less  sensitive  to  the  line  perturbations  that  occur  in  the  digitization  process. 
Furthermore,  the  processing  time  requirement  of  the  method  does  not  depend  on  the 
way  the  parameter  space  is  quantized.  In  section  2.2,  this  algorithm  was  applied  on  a 
sample  image  and  its  advantages  were  discussed. 


In  section  3,  an  extension  to  the  above  algorithm,  called  WLF  algorithm  was 
presented.  The  main  advantage  of  this  technique  is  the  ability  to  define  arbitrary  criteria 
for  search  for  lines.  This  was  shown  in  section  3.1  as  well  as  several  other  advantages. 
While  the  WLF  technique  needs  the  same  processing  time  as  the  Hough  algorithm, 
it  needs  less  memory  to  store  the  results.  The  processing  time  is  independent  of  the 
number  of  quantization  levels  in  the  parameter  space,  thus  for  higher  k  (the  number 
of  quantization  levels  for  0),  it  needs  less  time  than  the  Hough  algorithm. 


Further  experimentation  on  real  data  from  general  images,  i.e.  ones  with  nonbinary 
pixels,  is  needed  to  establish  the  processing  time,  storage,  and  accuracy  of  the  weighted 
line-finding  algorithm  as  a  road  —  and  general  mapping/ charting/geodesy  feature  — 
locator. 
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APPENDIX  A 

Taking  a  line  passing  through  a  pair  of  points  (Xj,  Vj )  and  (X2>  Y2),  its  equa¬ 
tion  will  be 


(a2)  Y  =  Y,  +  (Y,  -  Y2)  •  (X-X,)/(X,  -X2) 

And  the  slope  of  the  line  will  be: 

(a3)  m  =  (Yi  -Y2)/(Xl  -X2) 

Since  at  point  A,  Y  =  O,  therefore  the  distance  OA  will  be 

(a4)  OA  =  (X2Yj  -X,Y2)/(Y,  -Y2) 

Taking  the  normal  to  the  line  from  the  origin,  the  slope  of  the  normal  will  be 
(a5)  m  =-l/m  =  -(X,  -X2)/( Y,  -Y2) 

So  the  equation  of  the  normal  will  be 

(a6)  Y  =  mx  =  -(X1  -  X2)  •  X/(Y,  -  Y2) 

By  equating  the  (a2)  and  (a6),  we  can  find  the  coordinates  of  the  intersection  (X^Yj): 
(a7)  Xi  =  [  (Xj  Y2  -  X2Yj  )  •  (Y2-Y,)l  /  l  (X,  -X2)2  +  (Y,  -Y2)2  ] 

(a8)  Yj  =  [(X1Y2-X2Y1)  •  (Xi  -X2)  |  /[(Xi  -X2)2  +  (Y,  ~Y2)2  1 
and  the  length  of  the  normal  to  line  1 ,  r, ,  would  be  found  from: 

(a9)  r,2  =  X.2  +  Y2 

l  i  i 


( a  1 0)  r  j  =  |  Xj  Y2  -  X2 Y,  |  V(X,  -  X2)2  +  (Y x  -  Y 2)r/  [(X,  -  X2)2  +(YrY2)2l 
=  I  X,Y2  -XjY,  |V(X,  -X2)4  +  (Yj  -Y2)r 
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TABLE  1.  Accumulator  Array  of  Fig.  3  for  the  Hough  algorithm 
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TABLE  2.  Accumulator  Array  of  Fig.  3  for  the  line-finding  algorithm 


TABLE  3.  Accumulator  Array  of  Fig.  3  for  the  weighted  line-finding  algorithm 


'  ~  ru 

ru 

<3 

.  »■»  ru  «■* 

n 

in 

!  ru  n  n 

.  ru  -o  m  n 

m 

n 

O' 

.  ru  «r  —  ru  -< 

o 

ru 

H 

.nrunnmo-«'0«9’ru  —  n 

'•rnnc,j«r'0-on<r*onnn 

ra 

o 

'H^-innrU'O'Omnnnru 

O' 

ru 

O' 

—  runrvDDco«r(MtM 

o 

O' 

^  ^  ^  ru  ru  ru  ru  «■*  ru 

■  ru  ru  ru  « 

^4 

iH 

CO 

—  •^ru^runiru'O 

i  ru  ru  a  ru 

— 

ru 

N 

^  ru 

:  n  —  ru-« 

ru 

n 

o 

1 

1 

1 

1 

1 

1 

’  ru  ru  ru  ru 

■a- 

in 

1 

1 

1 

1 

1 

n  *i 

ru 

in 

<? 

1 

1 

1 

1 

1 

1 

^  w  n  n 

ru 

-c 

n 

^  ru  ru 

.  ru  c  «s-  — 

n 

Is* 

ru 

«  ^  ^  ^ 

'  -c  n  c  f; 

ru 

03 

^  ^  ^  ru  cu  ^ 

1 

4 

5 

6 

n 

0- 

in 

i 

i 

i 

i 

1 

2 
1 

1 

.  ru  ru  ru  -o 

CO 

o 

11 

ru 

c 

h  h  ^  »■<  ru  ru  ru  ru 

'  ru  ru  n 

< 

ru 

o 

.  m  o  in  os 

O' 

o 

.*  ru 

a 

■ 

cu^o^run^m>orvCDO"^ 

.  t  i 

ru  ru  ru  ru 

ru 

n 

c 

18 


I 


i 


TABLE  4.  Accumulator  Array  of  Fig.  4  for  the  Hough  algorithm 


TABLE  5.  Accumulator  Array  of  Fig.  5  for  the  line-finding  algorithm 


